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We develop a general Hamiltonian treatment of spontaneous four-wave mixing in a microring 
resonator side-coupled to a channel waveguide. The effect of scattering losses in the ring is included, 
as well as parasitic nonlinear effects including self- and cross-phase modulation. A procedure for 
computing the output of such a system for arbitrary parameters and pump states is presented. For 
the limit of weak pumping an expression for the joint spectral intensity of generated photon pairs, 
as well as the singles-to-coincidences ratio, is derived. 
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I. INTRODUCTION 

Quantum states of light, including single photons and 
entangled photon pairs, are a critical resource in quantum 
information processing. While traditional techniques for 
generating entangled photon pairs have relied on bulk 
nonlinear optical materials and waveguides mm, a less 
cumbersome micro-scale implementation is needed for an 
integrated optics setting. Designing such “on-chip” inte¬ 
grated optical components is likely a necessary step in 
the physical realization of optical quantum information 
processing. Recent advances in the fabrication of reso¬ 
nant optical microstructures 011 have yielded promising 
chip-based candidates for entangled photon pair genera¬ 
tion. One such implementation involves a micron-scale 
optical ring resonator side-coupled to a channel waveg¬ 
uide. By pumping the ring via the side channel, the 
nonlinear optical response in the ring may give rise to en¬ 
tangled photon pairs generated by spontaneous paramet¬ 
ric down-conversion (SPDC) or spontaneous four-wave 
mixing (SFWM). By operating at wavelengths close to 
the ring resonances, dramatic enhancement of nonlin¬ 
ear conversion efficiencies has been demonstrated with 
peak pump powers of only a few mW m- Photon 
pair generation from both SPDC and SFWM in silicon 
ring resonators CM. classical four-wave mixing 000 
and optical parametric oscillation |14j in silicon nitride 
and doped silica glass rings, as well as frequency comb 
generation in silicon nitride and aluminum nitride rings 
fTCUTB] have recently been experimentally demonstrated. 
There has been some theoretical study of SFWM and 
SPDC in these systems for the case of a continuous wave 
pump mum, and perturbative calculations have been 
performed which demonstrate how controlling the pump 
pulse duration affects spectral correlations in the gen¬ 
erated entangled photon pairs but so far these 

have neglected the quantum effects of losses on the gen¬ 
erated pairs due to scattering of photons in the ring. 
Parasitic nonlinear effects such as self- and cross-phase 
modulation (SPM and XPM, respectively), which may 
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become important for strong pumps, as well as other 
strong-pumping corrections have not been fully included 
in these theoretical treatments of photon pair generation. 

In this paper we present a Hamiltonian treatment of 
third-order nonlinear optical processes in such a system, 
focusing on SFWM while accounting for SPM of the 
pump field, as well as XPM between the pump, signal 
and idler fields. We include scattering losses in the ring, 
and show that they may be modelled as an effective ex¬ 
tra channel into which ring modes can emit. Probabilities 
corresponding to single-photon detection events wherein 
a photon’s pair partner is lost to scattering, which are 
not directly calculable in the absence of a formal loss 
mechanism, can then be studied. 

In the next section we begin by laying out the Hamil¬ 
tonian for the full ring-channel system and bath. Heisen¬ 
berg equations of motion for each field in the channel 
and ring are developed in section |III[ and the notion of 
incoming and outgoing fields is formalized. In section |IV| 
a strategy is developed to solve these equations for ar¬ 
bitrary parameters and pump pulses. Within our treat¬ 
ment the outgoing and incoming field operators are re¬ 
lated in the frequency domain by spectral response func¬ 
tions, determined by solving equations of motion for the 
ring operators alone. A specific perturbative example for 
weak pumps is worked out in section [VJ the results of 
which we then use to calculate the joint spectral inten¬ 
sity of generated photon pairs, as well as the singles-to- 
coincidences ratio. 


II. MODEL HAMILTONIAN 

We begin with a Hamiltonian for the channel and ring 
system as shown in Fig. |T| We take the channel to ex¬ 
tend from z = —oo to z = +00 and to accommodate three 
fields of interest: one pump with reference frequency up, 
as well as a signal and idler field with reference frequen¬ 
cies us and uj, respectively. Each of these fields couple 
linearly to the ring, with the coupling assumed to take 
place at a single point in space, 2 = 0. We further as¬ 
sume the ring geometry and field frequencies have been 
chosen such that the signal, idler and pump frequencies 
lie within narrow, well separated resonances. To begin 
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FIG. 1. Ring and channel microstructure geometry with la¬ 
belled incoming, outgoing and ring modes. Upper “phantom 
channel” represents scattering loss mechanism in the ring. 


with we focus on an individual pair of signal-idler fields, 
and discuss the case of multiple sets of signal-idler fields 
in later sections. 

Extending the Hamiltonian considered earlier m, 
within our model the fully quantum-mechanical Hamil¬ 
tonian can be split into four parts, 


H — H c } lanne i Hring A Hcoupling T H^ath • (1) 

wherein ff c (, anne i refers to the free Hamiltonian for the 
channel fields, H r i ng to the ring, H coup u ng to the linear 
channel-ring coupling, and Hbath to the modes into which 
photons can be lost, as well as their couplings to the ring 
modes. We then have for the channel Hamiltonian 


H, 


channel 


=£ 


hu>i 


dzipj(z)'ip J (z) 


J L 

ihvj 


dz 


dz 


i>j{z) - h.c. 


•( 2 ) 


In this expression (and elsewhere in this work) the index 
J ranges over the mode labels {P, S',/}. The channel 
fields can be decomposed into constituent modes 

with lowering operators Aj(k) according to 

Mz) = J (3) 

where kj is the reference wavevector corresponding to 
frequency wj for the field ipj(z), and is subtracted in the 
exponential to lead to a 'tpj(z) that is slowly varying in 
space. If the integrals in © ranged from —oo to oo, the 
channel operators il>j(z) would satisfy 


[i/;j(z),ipj’(z')}= 0, 

= 6(z-z')5jj'. 


( 4 ) 


While the integrals in ([3]) do not span the whole range of 
k, both because of mode cut-offs and because the same 
transverse mode of the waveguide might be broken into 
different spectral ranges to identify different fields ipj{z), 
we assume we work with pulses sufficiently long that Q 
is a good approximation. The coefficients vj in H c h a nnei 
refer to the propagation speeds for each field J, account¬ 
ing for group velocity dispersion between the different 
channel fields. While our model neglects group velocity 
dispersion between modes within a single channel field, 
our treatment can easily be generalized to include arbi¬ 
trary dispersion. 

Whereas in our model each field in the the channel 
consists of a continuum of modes, we approximate the 
fields in the ring as isolated individual modes, which is 
valid provided the ring resonances are sufficiently narrow 
[251 . The ring Hamiltonian can then be written in terms 
of discrete mode annihilation operators bj , 

Hring — ^ ' bujjb^jbj + HjVL- (5) 

J 

The latter term in the above expression is responsible for 
all the nonlinearity in this system. Since the pump field 
intensity in the ring resonator will be much larger than 
that in the channel, we assume that H^l contains only 
ring operators. In this work we focus on effects arising 
from the third-order nonlinear susceptibility y^ 3 ), taking 

Hnl =(hXbpbpbgb\ + h.c.^ + Tirjb^pb^pbpbp 

+h( (b^pbsbp + b^pbjbp) . (6) 

The first term in this expression leads to photon pair 
generation: two pump photons may be converted into a 
signal and idler photon pair. The second represents SPM 
of the pump, while the latter two correspond to XPM be¬ 
tween the pump and the signal and idler modes. SPM of 
the signal and idler modes is neglected, as the number of 
photons in those modes will be small compared to that 
in the pump mode. The nonlinear coupling coefficients 
A ,rj and ( are in general not independent, for they arise 
from the same underlying optical nonlinearity, and can 
be calculated from the mode fields [24] in the ring. A sim¬ 
ple estimate gives h\ ss 3(ftwp) 2 xl 3 V(4eon 4 f2 I . inff ), with 
r) = A/2 and £ = 2 A, where n is the refractive index of the 
ring and firing the volume of the ring mode, and eo the 
permittivity of free space, assuming only one component 
of the nonlinear susceptibility \is accessed uniformly 
in the ring. For ease of exposition we concentrate on 
degenerate SFWM, though the methods presented here 
may easily be generalized to include multiple pump fields, 
as well as effects arising from other nonlinearities, includ¬ 
ing SPDC. Our formalism also readily encompasses clas¬ 
sical four-wave mixing if the input state in each mode is 
treated as a coherent state. 

We adopt a point-coupling model for the interaction 
between the channel and ring. This approximation is 
justified provided the spatial extent of the interaction 
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region is short compared to all relevant wavelengths [25]. 
Our ring-channel coupling Hamiltonian is then 

//coupling E + h.c.) , (7) 


in which we have introduced the ring-channel coupling 
coefficients 7 j. These constants can be related to the 
usual self- and cross-coupling coefficients [25] used to 
characterize microring resonators |26l . 

In an experimental setting, loss due to scattering may 
significantly affect the yield and statistics of the gener¬ 
ated photons, manifested as a nonzero probability to de¬ 
tect a lone signal or idler whose pair partner was lost 
BM- We formally account for these effects by the in¬ 
clusion of Hb a th , which would in principle accommodate 
modes into which ring photons can scatter. Since the lo¬ 
cal field intensity is strongest within the ring resonator, 
we assume that the effects of loss there dominate those 
in the channel. We therefore model the bath Hamilto¬ 
nian as a set of continua of modes into which each ring 
mode can emit. As demonstrated in Appendix [A] this 
loss mechanism is equivalent to introducing an effective 
“phantom channel”, identical to our physical channel, to 
which each ring mode linearly couples. We therefore take 

Hbath — Hphantom channel ~t“ Hphantom coupling (3) 

where 


Hphantom channel 

At , 


E 


j 


Tiloj / dz(p'j(z)(l>j(z) 


+ 


ihuj f , l d(j)j(z) t . dcj)j(z) 


dz 




and 


H, 


(9) 

phantom coupling — E (^j^j(O) +h.c.) . ( 10 ) 


In the above expressions we have introduced effective 
phantom channel fields 4>j, satisfying the same commu¬ 
tation relations as Q. The effective velocities uj are 
allowed to differ from the vj, and the reference frequen¬ 
cies ujj are identical to those for the ring modes. In this 
manner a photon lost to scattering can be thought of as 
existing and propagating away in the effective loss chan¬ 
nel. When calculating physical quantities with our model 
we need only average over all such loss events. 


III. EQUATIONS OF MOTION 

Having constructed our Hamiltonian, we now adopt 
a Heisenberg picture and compute the equations of mo¬ 
tion for each operator. We then simplify the coupling 
description by introducing slowly-varying incoming and 
outgoing field operators, reducing the entire problem to 
solving a set of coupled nonlinear ordinary differential 
equations for the ring operators alone. 


A. Heisenberg equations 


Let to be the initial time. For each Schroedinger 
operator O = 0(to) we introduce the corresponding 
Heisenberg-picture operator 0(t) via 

0(t) = Ui(t,to)0(to)U(t,to), 

( 11 ) 

where the time evolution operator t/(t, to) 
{/(to, to) = 1 and 

satisfies 

ih ^ ^ -HUit, t 0 ). 

( 12 ) 

The equation of motion for O(t) is then 


^ = i- 

(13) 


We now proceed to calculate the equations of motion for 
each of the operators that appear in the Hamiltonian ([l]) . 
For the channel fields we obtain 

f d d \ 

\dt, +VJ lf' +IUJJ ) = ( 14 ) 

essentially a propagation equation with a driving term 
at z = 0. An analogous equation holds for the phantom 
channel fields (f>j. The ring pump operator is found to 
satisfy 

b P (t) (15) 

= -ryptMM) - */rp</>p(0,t) - 2i\*tf P (t)b s (t)b I (t). 

For the signal and idler modes in the ring, we obtain for 
the Heisenberg operators 

( + iu} S + iCb f p (t)b P {t)^bs(t) = (16a) 

- *7sV>s(0, t) - iv* s <Ps(0, t ) - iXbp{t)b P {i)b\{t), 

bi(t) = (16b) 

- *7^/(0, t) - 0, t) - iXbp(t)b P (t)bl(t). 


( ^ + + i(b\,{t)b P (t)j 


^ + iojp + 2irjb 1< p (t)bp(t) S j 


B. Incoming and outgoing fields 


The equations of motion (14) for the channel fields can 


easily be formally solved, giving 
= T j{z,t ) 

bj {t - */ VJ ) [9iz) - 0(z - vjt )] (17) 

vj 

where j(z, t) is any solution to the homogeneous version 
of ( fl4| , which must take the form ^jiz,t) = e~ luljt fiz — 
vjt), with / being any operator-valued function of one 
variable. We may substitute this expression evaluated 
at z = 0 into our ring operator equations (|15|) and (16 1 . 


Note that as a consequence of our point-coupling model, 
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this solution for i/jj(z,t) is discontinuous across z = 0 , 
a feature that survives even in classical applications of 
the model. It is dealt with 123 ESI by introducing the 
one-sided limits 

ipj(0 ± ,t)= lim ipj(z,t) (18) 

z->0± 

and defining ■!/> j(z,t ) at z = 0 via 

V>j(M) -t ^ + . (19) 

Using these definitions, we can relate each field at the 
coupling point z = 0 to the corresponding ring opera¬ 
tor and the field immediately left of £ = 0 (that is, the 
“incoming” field), 

Vv(o,t) -t - \ l — bj(t). ( 20 ) 

2 vj 

Each field operator immediately to the right of the cou¬ 
pling point (the “outgoing” field) can now be put in terms 
of the incoming field and the corresponding ring opera¬ 
tor: 


To avoid explicitly dealing with field discontinuities at 
the coupling point, we construct formal entities that we 
identify with incoming and outgoing fields. Define the 
incoming field ij)j< via 

ipj<(z,t) = ipj(z,t) for 2 < 0 (26) 

and extend it to z > 0 by demanding everywhere that 

( Q Q \ 

Q t + v J-g z +iu j Jip j< {z,t) = 0. (27) 

That is, we give to ipj < a false future to the right of the 
coupling point, which corresponds to the free evolution 
in the absence of any coupling to the ring. We also define 
the outgoing field t/jj > via 

i/jj > (z,t) = ipj(z,t) for z > 0 , (28) 

and similarly give it a false past by enforcing everywhere 

(J^ +VJ § Z + * WJ ) Vu>(z,*) = 0- (29) 


Vu(0 + ,f) = ipj(0 ,t) - —bj(t). (21) 

vj 


Using (20 1 (and the analogous expression for the phan¬ 
tom channel fields <j>j) in our ring operator equations (15) 
and ( 161 , we obtain 


+ Tp + iuip + 2ir]b^ P (t)bp(t)^j b P (t ) 

= - i'Yp^pfi-,t) - in*p(j)p(Q~ ,t) 

-2i\%(t)b s (t)b I (t) (22) 


for the pump, and 

^ ^+F s+ivs + i(b f p (t)b P (t)'jbs(t) (23a) 

= -*7s^s(0“,t) - i/4<As(0 _ ,t) - iXbp(t)b P (t)b\(t), 

^ ^ +T/ + <^p(t) 6 p(t)^fe/(t) (23b) 

= ,t) -f/x^s( 0 _ ,t) - i\bp(t)b P (t)bl(t) 


for the signal and idler. By formally solving for the chan¬ 
nel fields and self-consistently substituting the result into 
the equations of motion for the ring operators, effective 
damping terms Tj therein are introduced, 


where 


' j=Tj + Mj , 

(24) 

r ItjI 2 

rj = 2vj • 

(25a) 

2 uj 

(25b) 


The Tj damping constants each contain a contribution 
from the coupling to the physical channel, as well as a 
contribution from loss to the phantom channel. 


Since we will be interested in the properties of the gen¬ 
erated photons, which exit the ring and propagate to the 
right, all physical quantities of interest (such as photon 
pair generation probabilities) will involve the outgoing 
channel fields ipj>. Indeed, since these fields follow the 
equation of motion corresponding to free evolution, the 
field at large positive 2 (where measurements on the gen¬ 
erated light would occur) is completely determined by 
the behaviour at z = 0 : 

= e~ iujZ t Vj ip J> (Q, t — -/vj) for * > 0. (30) 


Our goal is therefore to construct a solution for these 
outgoing fields in terms of the incoming fields ij)j <, with 
which our initial pump state will be expressed. It suffices 
to do so at the origin, since the outgoing fields for z > 0 
are trivially related to those at z = 0 . 

Another advantage of the trivial free evolution dynam¬ 
ics obeyed by if>j< and Vu> is that each field can be 
Fourier decomposed in a simple way: 

j -^=aAn)e i ^e- i ^+^ t , (31a) 

1>j>(z,t)= j -^c J (K)e to e- i ^+^)‘, (31b) 

wherein the aj(n) and cj(k) are time-independent. 
These amplitudes are essentially the annihilation opera¬ 
tors for the modes that make up the i/jj< and ^j> fields 
on the respective half lines z < 0 and z > 0 , and obey 
the correct bosonic commutation relations 



S(k - k')Sjp, 


cj(k ), 



S(k - k')Sjj>. 


(32) 


Within our treatment calculations are most naturally 
done with respect to these annihilation operators in place 
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of the Aj(k) defined in equation ([3]). We may perform 
the exact same procedure for the effective phantom chan¬ 
nel fields, introducing 

<t>j<(z,t)= J (33a) 

<f>j>(z,t)= J -^=gj( K )e i **e- i < u '+ miJ *. (33b) 

The outgoing phantom channel fields (f>j > will contain 
all the photons lost to scattering or absorption. The in¬ 
coming phantom channel fields will also be important, as 


they will amount to an extra source of vacuum fluctua¬ 
tions injected into the ring, which are critical for seeding 
the spontaneous nonlinear processes that give rise to gen¬ 
erated signal and idler photons. 

It will be convenient to remove the rapidly varying 
phase from the incoming and outgoing fields, as well as 
the ring operators, defining the full envelope quantities 

Oj(t) = e iuijt Oj(t) (34) 

for each operator Oj(t). In terms of these quantities, our 
ring operator equations become 


(j^ + T P + 2ir]bp(t)bp(t) \ b P (t) = -*7 P V’p<(0,£) - i/r P ^ P< (0,f) - 2iX*bp(t)b s (t)bi(t)e~ zAt , (35a) 

(J^+rs+i$p(t)bp(t)j b s (t ) = -*7s'0 s< (O,t) - i/4^ s< (0,f) - iXb P (t)b P {t)b\(t)e lAt , (35b) 

+T/ + i(P P (t)b P (t)^J bi(t) = ^*7/'0/ < (O, t) - i/ij0j<(O,t) - iXb P {t)l P {t)b\(t)e lAt , (35c) 


where we have introduced the detuning parameter 

A = ws + w/ — 2uj p . (36) 

These ring equations, together with the transformation 
between the incoming and outgoing channel modes, 

VU>(M) = ip J< {0,t)- — bj(t), (37) 

vj 

constitute the fundamental equations of interest for the 
ring-channel system. 


IV. SOLVING THE EQUATIONS OF MOTION 


In the previous section, the problem of calculating the 
outgoing fields in terms of the incoming fields was re¬ 
duced to solving a set (35) of coupled ordinary differ¬ 


ential operator equations. These equations fully retain 
the quantum mechanical nature of all light modes, and 
may therefore be used to treat arbitrary injected pump 
states. In many experiments, however, the pump is a co¬ 
herent laser beam or pulse train well described by 

a classical function of time. As a first example we study 
the semi-classical versions of equations (35), in which the 


pump mode b P (t) is replaced by its expectation value: 


b P (t) -A P P (t) = {b P (t)). 


(38) 


At this stage we also implement the undepleted pump ap¬ 
proximation. In the ring pump equation (35a) the term 


multiplying bs(t)bi(t) accounts for the effect on the pump 
mode when a photon pair is generated. Neglecting such 
effects, we drop this nonlinear term in the pump equa¬ 
tion. We note, however, that pump losses are still ac¬ 
counted for within this approximation, as evidenced by 


r 


the coupling of the pump operator to the correspond¬ 
ing phantom channel field. “Undepleted” in this context 
therefore refers only to the nonlinear interaction; the ef¬ 
fects of linear coupling and loss are included in our de¬ 
scription of the pump. 

With the semi-classical substitution, the (undepleted) 
pump equation (35a) becomes an entirely classical and 


self-contained ordinary differential equation. In the ab¬ 
sence of SPM its solution is trivial, but in general it must 
be solved numerically. The signal and idler equations 
(f3o]>c) may be solved via a Green function strategy: 
these equations may be written in the form 


d_ 

dt 


bs(t) 

b\(t) 


= M(t ) 


bs(t) 

b\(t) 


D(t), 


where we have introduced the matrix 


M(t) = 


-T s -iC\P P (t)\ 2 -iX\p P {t)\ e' 


At 


iX* P P (t) 


0 —iAt 


-T I + iC\P P (t)\ 2 


as well as the driving term 

m = 


-n *si>s< ( 0 , t) - Vs^s< (°, t) 
( 0 , t) + (0 ,t) 


(39) 


(40) 


(41) 


Provided a 2 x 2 matrix G(t,t') may be found which 
satisfies 


dt 




for t > t! subject to initial conditions 
G(t\t')= [J J 


(42) 


( 43 ) 
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then the full solution to (39) is simply 


b s(t) 

b\(t) 


bs(t 0 ) 


b\(to) 


— G(t, to) 

+ f dt'6{t-t')G(t,t')D(t'). 


(44) 


The signal and idler ring operator dynamics are thus eas¬ 
ily transformed into a single 2x2 matrix ordinary differ¬ 
ential equation (421. Since the matrix G(t , t') is complex¬ 


valued (as opposed to operator-valued), this system is 
easily solvable numerically. The full effect of XPM and 
SPM, as well as strong pumps and nontrivial pump pulse 
envelopes can then be investigated. We intend to carry 
out such a general analysis in future communications. 


V. PERTURBATIVE SOLUTION 


While (42) must ultimately be integrated numerically, 


as an analytic example we now work out the weakly 
driven case, where we expect a perturbative approach 
to adequately capture the physics of photon pair gener¬ 
ation; we neglect XPM and SPM, setting r) = £ = 0. 
For simplicity we also take the loss-channel propagation 
speeds to equal those in the physical channel, uj = vj 
for each J, though our results (59 63) do not depend on 


this assumption; they depend only on the full quantities 
r j and Mj fl24l[25 1. 


In solving our fundamental equations the physics will 
be most intuitive working in the frequency domain. We 
define for each ring operator the corresponding Fourier 
amplitudes: 


so that 


bj(t) = J -^-bj(K)e~ iKVjt (45) 

bj{K) = vj J JLbjWe™*. (46) 


A. Channel transformations 

Consider first the (semi-classical and undepleted) 
pump equation (35a). Let our initial time be to- The ho¬ 


mogeneous contribution to to the solution of this driven, 
damped, linear ordinary differential equation can be 
dropped by assuming there is no field in the ring at to- 
Remaining then is the particular solution /3p(t), which 
in the frequency domain will satisfy 


—IKVp + 1 p 


(47) 


where we have taken (dp(n)) = 0, since we assume there 
is no incoming field in the phantom channel. This solu¬ 
tion is exact within the semi-classical undepleted pump 
approximation. 


We may perform the same procedure to find bs(n) and 
6 /(k) to zeroth order in the nonlinear coupling. In par¬ 
ticular, for the idler we obtain to zeroth order 


t(o)/ , -*7ja/(K) - ip^d^n) 

b i ( K ) =-:-—p-■ 

—IKVl +I/ 


(48) 


Inserting this and the pump solution into the frequency- 
domain version of our signal equation ( |35b[ ) , we find that 
to first order in the nonlinear coupling &s(«) satisfies 

(~iKv s +r s )bs(K) = —i'Ysas(n) - ip,* s ds(n) 

—i\vs [ dn'F s(k,k') b\\n') ,(49) 


where 


F s (k,k) = 


— J dni J dK 2 /3p{Ki)/3p(K2) 

xS^KVs — K\Vp — K2 Vp + k'vi + A) 


(50) 


Solving for bs( k) and substituting the result into the 
frequency-domain version of (37), 


c s (k) = a s (n ) - —&s(k), 
v s 


(51) 


we obtain for the outgoing signal field amplitudes in the 
channel 


cs(k) = J dn' 


Qss{k, n')as(n') + Pss(k, n^dsin 1 ) 


+<?S/(k, K')a\(K') +Psi(k, n')d\(n') 


, (52) 


where we have introduced the spectral response functions 


, -invs -T s + M s , 

qss{K, K)= -:-=- d(K-K), (53a) 


-invs + Tg 
t m ~lSP*shs ,, A 

PSS{K, K )= — - =-0{K - K ), 


9s/(k,« , )= 

Ps/(k,k , )= 


-invs + Tg 

-i'ys'yAF s (K,n') 
{-invs + Ts){in'vi + T/) ’ 
-ilsPiXFsjn, n') 

( inv s + r s ){iK , v I + r 7 ) 


(53b) 

(53c) 

(53d) 


Using an identical procedure, we may write down the 
channel transformation for the idler modes. We obtain 


Cj(/S) = Jd 


qn(K, K')a/(«/) +ph(k, k')^/(«') 


+qis{n, K ') a s( K ') + Pis{k, K^dlin') 


'Wt t’\ 


,(54) 


where the quantities qi x and pi x for x = {S', /}, as well 
as Fi, can be obtained by swapping the labels I in 
(53) and (50). Analogous functions can also be derived 
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for the transformation between the incoming fields and 
the outgoing loss-channel fields. 

From these expressions emerges the outline of a phys¬ 
ical picture of the system’s dynamics in the regime of 
weak nonlinearity. Within this perturbative approxima¬ 


tion, equations (521 and (53) describe the transformation 


between the incoming and outgoing signal field modes in 
the channel. This formalism is reminiscent of the conven¬ 
tional treatment of an optical beamsplitter, in which the 
outgoing and incoming mode annihilation operators are 
related by a simple linear transformation C2ZI- For our 
system there are a continuum of incoming and outgo¬ 
ing modes in the channel, with the ring playing the role 
of a “nonlinear beamsplitter” that scatters pump pho¬ 
tons into the signal and idler modes. As evident from 
the Dirac delta distribution in (|53a|) , the outgoing signal 


modes are correlated only with incoming signal modes of 
equal energy. On the other hand, from (53c) we see that 


the outgoing signal modes are correlated with incoming 
idler modes having a range of energies determined by the 
pump-dependent function Eg (re, re'), which describes an 
average over all energy-conserving ways for two pump 
photons to be converted into a signal and idler photon 
pair. 

It is important to note that while this perturbative 
solution treated the pump semi-classically, it is straight¬ 
forward to retain the pump’s quantum-mechanical na¬ 
ture. Keeping the pump modes described quantum- 
mechanically allows us to treat, for example, the case of 
a two-photon Fock state in the pump producing a signal 
and idler pair. More generally, this solution can read¬ 
ily encompass pump states with zero expectation value 
for the electric field. Our semi-classical case refers to 
the more common laboratory setting, involving a coher¬ 
ent pump laser pulse, which may be treated by taking 
the initial state of system to be a coherent state in each 
pump mode, and vacuum elsewhere: 


\^(to)} = e 


_/ dKQp(/t)a|,(K)-h.c. 


vac), 


(55) 


where ap(n) is the amplitude in mode re of the incoming 
pump field. Within this state our semi-classical pump 
amplitude (3 P (re) is simply 


/?p(k) = 


-*7 p 


—invp + Tp 


7=-ap(K). 


(56) 


Also evident from the structure of our perturbative so¬ 
lution is how the generalization to multiple signal and 
idler fields would play out. At the level of a weak unde¬ 
pleted pump, the presence of extra signal and idler modes 
would not affect the individual q xx > (re, re') and p xx > (re, re') 
functions. Each additional signal and idler mode would 
have its own set of response functions, formally identi¬ 
cal to (53), with the appropriate coupling constants and 


B. Physical quantities 

Having derived the transformations between the in¬ 
coming and outgoing channel field operators in the per¬ 
turbative regime, we are ready to calculate physical quan¬ 
tities. We will represent the initial state in the form 


of equation (55) as a classical, resonant Gaussian pump 


pulse of duration a, 


crp(re) = a P e 


= tv„ p -( ku W 2 ) 2 


(57) 


with vacuum in all other modes. 

Our focus will be on those physical quantities mea¬ 
surable with frequency-resolving photodetectors located 
at large positive z in the outgoing channel. Photode¬ 
tection probabilities are then given by the appropriate 
Glauber formula involving the expectation value of prod¬ 
ucts of outgoing channel field creation/annihilation op¬ 
erators [27]. The probability of detecting a signal and 
idler photon pair with wavevectors offset from the refer¬ 
ence wavevectors res and re/ by re and re', respectively, is 
proportional to 

$(re,re') = (c^(re)c|(re , )cs(re)c/(re')), (58) 

where the expectation value is taken with the state |vac), 
the vacuum of the incoming channel field operators. This 
is exactly the joint spectral intensity of the signal and 
idler fields. Using our perturbative response functions 
(53) we obtain (to leading order in A, in keeping with the 
order of the perturbative solution) 


<I>(re, k') = 


|A| 2 | 7 5| 2 |7/| 2 |*M«X)| 2 


((reu s ) 2 


r|)((re'u /) 2 + r 2 )' 


(59) 


The shape of this function is determined by the pump 
through Fs(k,k') - in particular, the pump pulse band¬ 
width determines the width of this two-photon distribu¬ 
tion, as well as the degree of spectral correlation between 
the two photons in the pair. 

Also evident from this expression is the effect of loss 
on <F(re, re'). Equation (59) is in formal agreement with 
earlier work m where loss was neglected, but with Tj 
now replaced by Tj. That is, the inclusion of loss is 
completely encapsulated in the total effective damping 
rates Tg and T/, which affect the widths of the overall 
Lorentzian envelopes in the denominator, as well as the 
transfer of pump energy from the channel to the ring as 
described by equation (561. While the phantom channel 


fields do contribute to the vacuum fluctuations that seed 
the pair generation process, as evidenced by the presence 
of the p xx i functions in (52), the overall yield of pairs 


linewidths. 


that exit to the physical channel only depends on loss 
insofar as the ring linewidths are affected. Plots of the 
normalized joint spectral intensities for a 100 ps pulse 
incident on a ring without and with loss (Tj = 10 GHz 
and Tj = 20 GHz, respectively) in the regime where the 
pulse duration is comparable to the lossless ring photon 
lifetimes are illustrated in Fig. <§. In this regime, in the 
absence of loss, the signal and idler energies are nearly 
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FIG. 2. Normalized joint spectral intensity function for a 100 
ps pulse incident on a ring (a) without loss, Tj = 10 GHz 
and (b) with loss, Tj = 20 GHz. Propagation speeds for all 
modes are taken to be vj = 15 cm/ns. 


FIG. 3. Normalized joint spectral intensity function for a 1 
ns pulse incident on a ring (a) without loss, Tj = 10 GHz 
and (b) with loss, Tj = 20 GHz. Propagation speeds for all 
modes are taken to be vj = 15 cm/ns. 


uncorrelated, as evidenced by the circular shape of the 
joint spectral intensity. The inclusion of the phantom 
channel with coupling strengths equal to those for the 
physical channel doubles the total effective damping rate, 
leading to significant broadening of the overall distribu¬ 
tion, and introduces some degree of spectral correlation. 
Fig. § illustrates the opposite regime, in which a longer 
1 ns pulse is employed to yield entanglement between the 
signal and idler photons. Loss similarly leads to broaden¬ 
ing of the overall signal and idler linewidths, but does not 
alter the degree of spectral correlation, indicated by the 
breadth of the distribution in the direction orthogonal to 
the antidiagonal. 

Though loss can be fully accounted for in the joint 
spectral intensity as corrections to the ring linewidths, 
the inclusion of the phantom channel introduces impor¬ 
tant quantum effects in the statistics of the outgoing pho¬ 
tons in the physical channel. These arise from the prob¬ 
ability for one photon in a generated pair to be lost to 


scattering in the ring, with only its partner subsequently 
available for detection. The total rate of pair detection 
events - the “coincidences rate” - is simply proportional 
to 

P'coincidences = J d/t J dhv ^(ft, K ). (60) 

This quantity must be distinguished from the total sin¬ 
gles probability, 

(4( K )3^(« , )C5(K)5/(K , )) (61) 

+(c\( k )9s( k ') c i( k )9s(k')) , 

which represents the detection of a signal or idler photon 
whose entangled partner was lost to scattering in the ring. 
Of particular interest is the ratio between this quantity 


P> 


singles 


= I dn I dn' 
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(the “singles rate”) and the coincidences rate, 


r = 


P f 


singles 


Pcoincidences 


(62) 


The inclusion of a loss mechanism in our treatment allows 
us to calculate this quantity, which in previous theoret¬ 
ical treatments would always equal zero. A more realis¬ 
tic estimate of r is especially important to applications 
wherein one of the photons in a pair is used to trigger 
another branch of the experiment that requires the trig¬ 
ger photon’s pair partner. Actuating a detector off of a 
photon whose partner was lost would typically be unde¬ 
sirable, making it necessary to understand the origin and 
magnitude of this effect. In the perturbative regime, we 
find 


TgMj + T/Mg 

r s r 7 


(63) 


intend to present further analyses of this regime in future 
works. 

Appendix A: An equivalent channel model for loss 

To model scattering losses in the ring, we imagine cou¬ 
pling each ring mode J to a continuum of bath modes 
described by annihilation operators Dj(k). The bath 
Hamiltonian, including the coupling to the ring modes, 
would be of the form 

H-bath £ dnhuj j(k)D*j(k)D j(k) 

j 

+ ^b\ J dnhp,j(K)Dj(K) + h.c^j (Al) 

where ujj(k) are the bath mode frequencies, and fij(n) 
are coupling coefficients to the corresponding ring modes. 
We now introduce a fictitious position parameter z, al¬ 
lowing us to define field operators 4>j(z) in the usual way, 


This expression is entirely independent of the pump and 
nonlinear coupling strength. In actual experimental im¬ 
plementations, the ring-channel system is often operated 
at or near the critical coupling regime, where the lin¬ 
ear transmissivity through the device is minimized eh. 
Within our treatment this occurs when the coupling rates 
between the ring and the phantom channel equal those 
between the ring and the physical channel: Tj = Mj for 
each J. In this regime (631 gives r = 2. Other losses 
present in an actual experiment can be expected only to 
increase this value, so r = 2 can be taken as the theoret¬ 
ical minimum for lossy, critically coupled ring resonator 
systems. 


VI. CONCLUSION 


We have developed a general Hamiltonian treatment 
of spontaneous four-wave mixing in a microring res¬ 
onator side-coupled to a channel waveguide. The inclu¬ 
sion of losses in the ring permitted the calculation of the 
singles-to-coincidences ratio r, for which a lower bound 
of r = 2 was derived in the critically coupled weak pump¬ 
ing regime in the absence of other losses. In this regime 
the inclusion of loss was found to affect the joint spectral 
intensity distribution of generated photon pairs strictly 
through corrections to the linewidths of the ring reso¬ 
nances. As such, when the pump pulse duration is com¬ 
parable to the lossless ring photon lifetimes, the addition 
of loss introduces significant broadening and some spec¬ 
tral correlation to the generated photon pairs. For longer 
pump pulses loss merely broadens the overall joint spec¬ 
tral intensity, but does not appreciably change the de¬ 
gree of entanglement. In the strongly pumped regime, 
wherein self- and cross-phase modulation may become 
important, a Green function procedure for numerically 
deriving relevant physical quantities was presented. We 


M*) = j dne i ^ K ~ Kj ^ z Dj(k), (A2) 

where kj is the reference wavevector for the mode J. 
Provided the function jxj{n) is sufficiently smooth and 
slowly varying over the relevant range of k, we can take 

J dKhfijWe**-*''** » hnj5{z), (A3) 

where /xj is a constant. The coupling term in Hb at h then 
becomes 

J dKhflj(K)Dj(K) + h.c^j 

~ ^2 0) + h.c.'j . (A4) 

j 

If we now expand ujj(k) to first order in k about the 
the appropriate reference frequency, 

Wj(fi) = uj(kj) + (k - Kj)uj, (A5) 

where uj is the effective propagation speed in the phan¬ 
tom channel, 

uj= (wly (A6) 

then the full bath Hamiltonian becomes 


Hbath — ddphantom channel A Hphantom coupling i (A7) 

fined in ([9]) and (101. 


ACKNOWLEDGMENTS 

We acknowledge support from the Natural Sciences 
and Engineering Research Council of Canada. Thanks 
are also due to Nicolas Quesada for helpful discussions. 






10 


[1] H. Takesue and K. Inoue, Phys. Rev. A 70, 031802(R) 
(2004). 

[2] M. Fiorentino, S. M. Spillane, R. G. Beausoleil, T. D. 
Roberts, P. Battle, and M. W. Munro, Opt. Express 15, 
7479 (2007). 

[3] Q. Xu, D. Fattal, and R. G. Beausoleil, Opt. Express 
16, 4309 (2008). 

[4] F. Xia, M. Rooks, L. Sekaric, and Y. Vlasov, Opt. Ex¬ 
press 15, 11934 (2007). 

[5] M. Ferrera, L. Razzari, D. Duchesne, R. Morandotti, 
Z. Yang, M. Liscidini, J. E. Sipe, S. Chu, B. E. Little, 
and D. J. Moss, Nature Photonics 2, 737 (2008). 

[6] M. Ferrera, D. Duchesne, L. Razzari, M. Peccianti, 
R. Morandotti, P. Cheben, S. Janz, D.-X. Xu, E. Lit¬ 
tle, S. Chu, and D. J. Moss, Opt. Express 17, 14098 
(2009). 

[7] A. C. Turner, M. A. Foster, A. L. Gaeta, and M. Lipson, 
Opt. Express 16, 4881 (2008). 

[8] S. Azzini, D. Grassani, M. J. Strain, M. Sorel, L. G. Hclt, 

J. E. Sipe, M. Liscidini, M. Galli, and D. Bajoni, Opt. 
Express 20, 23100 (2012). 

[9] S. Azzini, D. Grassani, M. Galli, L. C. Andreani, 
M. Sorel, M. J. Strain, L. G. Helt, .]. E. Sipe, M. Lis- 
cidini, and D. Bajoni, Opt. Lett. 37, 3807 (2012). 

[10] S. Clemmen, K. Phan Huy, W. Bogaerts, R. G. Baets, 
P. Emplit, and S. Massar, Opt. Express 17, 16558 (2009). 

[11] E. Engin, D. Bonneau, C. M. Natarajan, A. S. Clark, 
M. G. Tanner, R. H. Hadfield, S. N. Dorenbos, V. Zwiller, 

K. Ohira, N. Suzuki, H. Yoshida, N. Iezuka, M. Ezak, 
J. L. O’Brien, and M. G. Thompson, Opt. Express 21, 
27826 (2013). 

[12] D. Grassani, S. Azzini, M. Liscidini, M. Galli, M. .]. 


Strain, M. Sorel, J. E. Sipe, and D. Bajoni, Optica 2, 88 
(2015). 

[13] R. Wakabayashi, M. Fujiwara, K. Yoshino, Y. Nambu, 
M. Sasaki, and T. Aoki, Opt. Express 23, 1103 (2015). 

[14] J. S. Levy, A. Gondarenko, M. A. Foster, A. C. Turner- 
Foster, A. L. Gaeta, and M. Lipson, Nature Photonics 
4, 37 (2010). 

[15] Y. Okawachi, K. Saha, J. S. Levy, Y. H. Wen, M. Lipson, 
and A. L. Gaeta, Opt. Lett 36, 3398 (2011). 

[16] H. Jung, C. Xiong, K. Y. Fong, X. Zhang, and H. X. 
Tang, Opt. Lett. 38, 2810 (2013). 

[17] J. Chen, Z. H. Levine, J. Fan, and A. L. Migdall, Opt. 
Express 19, 1470 (2011). 

[18] R. M. Camacho, Opt. Express 20, 21977 (2012). 

[19] L. G. Helt, Z. Yang, M. Liscidini, and J. E. Sipe, Opt 
Lett. 35, 3006 (2010). 

[20] Z. Yang, P. Chak, A. D. Bristow, H. M. van Driel, R. Iyer, 
J. S. Aitchison, A. L. Smirl, and J. E. Sipe, Opt. Lett. 
32, 826 (2007). 

[21] Z. Yang and .]. E. Sipe, Opt. Lett. 32, 3296 (2007). 

[22] L. G. Helt, J. E. Sipe, and M. Liscidini, Opt. Lett. 37, 
4431 (2012). 

[23] M. Liscidini, L. G. Helt, and J. E. Sipe, Phys. Rev. A 
85, 013833 (2012). 

[24] J. E. Sipe, N. A. R. Bhat, P. Chak, and S. Pereira, Phys. 
Rev. E 69, 016604 (2004). 

[25] J. Heebner, R. Grover, and T. Ibrahim, Optical Microres¬ 
onators: Theory, Fabrication and Applications (Springer, 
2008). 

[26] P. Chak, S. Pereira, and J. E. Sipe, Phys. Rev. B 73, 
035105 (2006). 

[27] C. Gerry and P. Knight, Introductory Quantum Optics 
(Cambridge University Press, 2004). 



